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Thermal balance of the jet in the source SS433 is considered with account of radiative 
and adiabatic cooling, and different heating mechanisms. We consider jet heating by the 
inverse Compton effect of coronal hard X-ray quanta on jet electrons, the influence of 
shock wave propagation along the jet, and jet kinetic energy transformation into heat via 
Coulomb collisions of jet and corona protons. The most important heating mechanism for 
the source SS433 turns out to be Coulomb collisions of jet particles with the surrounding 
medium. 
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1. Introduction 

SS433 is a unique massive X-ray binary system with precessing relativistic jets. It is 
situated at a distance of approximately 5 kpc sa 1.5 x 10 22 cm, nearly in the galactic 
plane. The optical companion V1343 Aquilae was first identified i n th e survey of 
stars exhibiting H a (656 nm) emission by Stephenson and Sanduleak^^. It is one of 
the brightest stars in the Galaxy, the bolometric luminosity of the object assuming 
isotropic radiation is Lt> i ~ 10 40 erg/s^l. SS433 is a close binary system with an 
orbital period of 13.1 days'" 1 . The uniqueness of this source comes from existence 
of narrow oppositely directed subrelativistic jets, which emit red and blue-shifted, 
periodically variable lines. A commonly accepted model of this object suggest that 
a continuous regime of supercritical accretion of gas onto the relativistic star is 
maintained. Thus, a supercritical accretion disk forms, together with narrow jets of 
gas propagating perpendicular to the disk plane from the central regions of the disk 
and having the relativistic speed of 0.26c. The optical star fills its critical Roche 
lobe, providing a powerful and almost continuous flow of gas into the region of the 
relativistic star at a rate of ~ lO~ 4 M0/yr^. 

The INTEGRAL observations of SS433 discovered X-ray emission in the range 
3 to 90 keV, which remains to be visible over the whole period of the binary. That 
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Fig. 1. Schematic picture of the source SS433 



leads to suggestion of the presence of an extended hot region (corona) in the central 
parts of the accretion disk^. Modeling of this spectrum by Monte-Carlo simulations 
provided an information about the physical parameters of this system. Here we 
analyze a thermal balance of the jet, and investigate sources of its heating, needed to 
maintain the temperature of the jet on the level, imposed by the X-ray and optical 
observations. The detailed description will be published inl 



2. Monte-Carlo modeling 

We accept, that this binary system consists of an optical star and a black hole, 
surrounded by an accretion disk, and a hot corona, with a couple of jets. We have 
considered the accretion disk within the radius 1.5 • 10 12 cm. Far from the central 
parts of the disk standard accretion disk theory ^ can be used to estimate its 
geometrical thickness h = r tan 9d ^ k , and the disk opening angle Odisk ~ 2 degrees. 
As follows from observations, the jets have a conical shape wi th t he opening angle 
0j e t ~ 1.2 degrees. All parameters are listed in the Fig Q] from US ro = 10 11 cm. 

The density profile n ~ r~ 2 , and the adiabatic temperature profile T ~ g 2 / 3 ~ 
r — V 3 5 have been accepted. It was obtained^Sthat th e best coincidence of the model 
with observational takes place at the adiabatic temperature distribution along the 
jet. In the figure Figj2] the best fit SS433 spectrum in the range from 3 to 90 keV 
is presented. The spectrum corresponds to precessional moment, when the angle 
between jet axis and the line of sight is equal 60 degrees and the disk is maximally 
face-on. 



3. Thermal balance of the jet 

The equation for the energy balance of the matter in the stationary jet with account 
of heating, and adiabatic and radiative cooling is written as^ 
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Fig. 2. The results of the simulations for different values of accretion disk effective temperature: 
6.0 ■ 10 5 K (solid line), 9.4 • 10 5 K (dash-dotted line), 1.1 ■ 10 6 K (dotted line). The best fit occurs 
when T B f = 6.0 ■ 10 5 K. r cor = 0.24, T cor = 19 keV, r cor = 6.4 ■ 10 11 cm, Mj Bt = 4 • 10 19 g/s. 



3 V jet k dT _ Vj et ^ 1 dp | 

2 Trip dr m p p dr 

Here e is the internal energy of the gas unit mass, q includes radiative cooling, 
and different types of heating, p is a pressure, V = 1/p is a specific gas volume. 
Bremstrahlung losses (erg/g/s) are determined by the expression, see^J 



*" = / ^ dVd6 = 16 (I) hr£r^^^ = M52.10-3n iet VT. (2) 

Solution of ([1} for q = qff, obtained in^, have shown a very rapid temperature drop. 
The value of T drops formally to zero at r as 1.5ro, what is in contradiction with 
observational data. In ^ all energy losses have been used in calculations, including 
bound-free, and bound-bound transitions from^. 

Jet heating by Compton interaction with hot photon gas emitted by the corona 
and the disk was investigated, and it was shown, that it is much less than the 
radiation losses ^ 



3.1. Heating by shock waves 

The shock wave generation may happen during the jet propagation through the 
surrounding gas. It was shown in ^ that due to the entropy increase in the shock, 
the average heating rate due to shock dissipation is determined as 



kT 



( 7 - l)m p M P 1 7 + 1 7+1 J +7ln U + l + M 2 ( 7 + 1)JJ (3) 
the value of At=l sec was taken equal to the characteristic time of SS433 variability 
in the X-rays^. The Mach number is decreasing due to transformation of the kinetic 
flow energy into heat in the shock.. The time averaged equation for the Mach number 
change along the jet is written as 
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1.4 ifi 



Fig. 3. The temperature profiles along the jet with account of heating by the shock waves, and 
full radiative losses. The dash line corresponds to M=2, dash-dot line corresponds to M=3, dash- 
dot-dot line corresponds to M=5. The dot line corresponds to the case without shock heating, and 
the full line corresponds to adiabatic expansion, tq = 10 sm. 
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The temperature distribution with account of shock heating and radiative losses is 
written as 



1 dT _ T dp a | m p / q sw q br \ 
7 - 1 dr po dr k \Mc s v jet J ' 

The solution of Q , ([5]) , for the shock formation at the disk origin, is shown in Figj3] 
It may be seen that the shock is rapidly dissipating. The shock heating could be 
important only if the shocks are originating in the jet all along its length^. 



3.2. Heating by Coulomb collisions 

Transformation of the jet kinetic energy into the heat happens when the protons 
from the corona and circumstellar gas are entering the jet. During collisions these 
protons are thermalized, transforming their relative kinetic energy into the heat, 
decreasing the kinetic energy of the jet. Maximal heating happens when the proton 
entering the jet loose its whole, and remains in the jet. In this case the rate of 
heating is determined as q C oii,max = v 2 et (v t h /2d jet )(n cor /n jet ) — r -1 . In reality, the 
proton may go through the jet, loosing only part of its energy. For the jet with 
above mentioned parameter coulomb collisions with protons of the jet may extract 
only small part of the energy because of small cross-section, at which the mean free 
path of the proton in the jet greatly exceeds its radius. Introducing t co u, equal to 
a ration of the proton energy lost in the jet, we have^ 
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Fig. 4. Temperature profile in the jet with account of the radiative cooling, and heating by the 
coulomb collisions, with a = (r/rrj) — 3 ^ 2 (dash line), a = (r/ro) — 1163 (dash-dot line). The full line 
corresponds to pure adiabatic cooling. 



2 f daAE If, , /Xd\ , s 

r coll = n jet d jet J , « 2-10- 4 , at , / ASda = 2^ In — . (6) 

" i P u jet " l P v jet/ z 'J \ Pp J 

Here p p = e 2 /m p v 2 et 2.1 • 1(T 15 , r = 6.4 • 10 11 , n jet = n jeti0 /41 2 • 10 13 cm- 3 , 
d jet = 2r tan (0 jet /2) 1.3 • 10 10 , In f ps 25. It follows that the jet heating due 

to coulomb collisions is much smaller than the radiative cooling ss 6 • 10 -3 . 

A presence even of a small magnetic field changes the situation. When a girora- 
dius Th of the proton is smaller than the radius of the disk the proton may spend 
larger time inside the jet and can loose a larger portion of its energy. The equality 
r L = Tj e t is reached at Bj et = '"^ CT ' h w 1.3 ■ 1CP 4 Gs, at the jet origin. In reality 
the magnetic field is probably much larger, so we may suggest that the proton is 
loosing a considerable part of its energy inside the jet. Here we try to find what 
part of the energy the proton should loose inside the inner, X-ray jet, to obtain the 
temperature distribution close to the adiabatic one. 

We choose the heating by collisions in the form q co ii(r) — q C oU.max(r) * a(r). 
Several examples for the temperature dependence at different choice a{r are given 
in FigHl The case with a = (r/ro) -1163 gives a good coincidence with the adiabatic 
curve up to 10ro, where the main source of the radiation is situated. It is clear that 
the whole adiabatic curve may be fitted with a more complicated function a = r). 

The optical jet in SS433 is situated at distances r ~ 10 14-15 cm (10 3_4 ro), and 
has almost constant temperature T ~ 10 4 ®. Because of radiative losses, and adia- 
batic expansion, a heating is necessary to support the constancy of the temperature. 
We consider the same mechanism of the optical jet heating. The protons come from 
the circumstellar gas with adiabatically decreasing temperature T ~ p 2 ^ 3 outside 
the outer radius of the corona r = r cor = 6,4 • 10 11 cm. Inside the corona the 
constant temperature is accepted To = 2.2 • 10 8 K. The density profile is taken the 
same inside and outside the corona n ~ r~ 2 , with the the density at the bottom 
of the corona n cor fi = 4.3 ■ 10 12 cm^ 3 . It was obtained in^l, that the temperature 
profile all over the jet may be reproduced by the choice of a(r), see Figj5]from^. 
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Fig. 5. Temperature profile in the jet with account of the radiative cooling, and heating by the 
coulomb collisions, with a(r) ~ r e3 ,r < r CO r\a(r) ~ r -0.48^ r y r C or (dash line); the full line 
corresponds to the analytic function a(r), which gives the best fit to the adiabatic curve in the 
X-ay jet, and constant temperature optical jet; dot line corresponds to a pure adiabatic expansion. 

4. Conclusions 

The most effective heating mechanism of the jet in SS 433 is connected with kinetic 
energy losses by collisions with surrounding matter (protons), in presence of a very 
moderate magnetic field. Collisions with protons from circumstellar gas may support 
T ~ 10 4 in the optical jet. Shocks distributed over the jet may give an input into 
heating. The losses of the jet kinetic energy to its heating are relatively very small 
1CP 4 , and the velocity change along the jet is hardly observable. 

Acknowledgments 

This work was partially supported by the RFBR grant 11-02-00602, the RAN pro- 
gramm P20, and RF President Grant NSh-3458.2010.2. 

References 

1. G.S. Bisnovatyi-Kogan, Stellar physics. Vol.1: Fundamental concepts and stellar equi- 
librium. Berlin, Springer (2001). 

2. G. S. Bisnovatyi-Kogan, Yu. M. Krivosheyev, Astron. Zh. (2011) (in press) 

3. A.M. Cherepashchuk, MNRAS, 194, 761 (1981). 

4. A.M. Cherepashchuk, A. A. Aslanov, V.G Kornilov, Sov. Astron., 26, 697 (1982). 

5. A.M. Cherepashchuk, R.A. Syunyaev, S.N. Fabrika, et al., Astron. Astrophys., 437, 
561 (2005). 

6. K. Davidson, R. McCray, ApJ, 241, 1082 (1980). 

7. S. Fabrika, ApSS Reviews, 12, 1 (2004). 

8. A.B. Kirienko Astronomy Letters, 19, 11 (1993). 

9. E.V. Koval', N.I. Shakura, in Proc. 23rd ESLAB Symp. on Two-Topics in X-Ray 
Astronomy, Bologna, Italy, 13-20 September, p.479 (1989), ESA SP-296. 

10. Yu. M. Krivosheyev, G. S. Bisnovatyi-Kogan, M. Cherepashchuk and K. A. Postnov 
MNRAS, 394, 1674 (2009). 

11. B. Margon, Annual review of astronomy and astrophysics. 22, 507 (1984). 

12. N.I. Shakura, R.A. Sunyaev, AAp 24 (1973) 337. 

13. C.B. Stephenson, N. Sanduleak, ApJ Suppl, 33, 459 (1977). 



